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Central limit theorem for supercritical binary homogeneous 
Crump-Mode-Jagers processes 

Benoit Henry^’^ 


Abstract 

We consider a supercritical general branching population where the lifetimes of individuals are 
i.i.d. with arbitrary distribution and each individual gives birth to new individuals at Poisson times 
independently from each others. The population counting process of such population is a known 
as binary homogeneous Crump-Jargers-Mode process. It is known that such processes converges 
almost surely when correctly renormalized. In this paper, we study the error of this convergence. 
To this end, we use classical renewal theory and recent works [Iliiiig on this model to obtain the 
moments of the error. Then, we can precisely study the asymptotic behaviour of these moments 
thanks to Levy processes theory. These results in conjunction with a new decomposition of the 
splitting trees allow us to obtain a central limit theorem. 
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1 Introduction 

In this work, we consider a general branching population where individuals live and reproduce inde¬ 
pendently from each other. Their lifetimes follow an arbitrary distribution Py and the births occur 
at Poisson times with constant rate b. The genealogical tree induced by this population is called a 
splitting tree [miiniiii] and is of main importance in the study of the model. 

The population counting process Nt (giving the number of living individuals at time t) is a 
binary homogeneous Crump-Mode-Jagers (CMJ) process. Crump-Mode-Jagers processes are very 
general branching processes. Such processes are known to have many applications. For instance, in 
biology, they have recently been used to model spreading diseases (see [2010]). Another example of 
application appears in queuing theory (see [l8] and [E]). 

In |19] . Nerman shows very general conditions for the almost sure convergence of general CMJ 
processes. In the supercritical case, it is known that the quantity e~^^Nt, where a is the Malthusian 
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parameter of the population, converges almost surely. This result has been proved in [21] using 
Jagers-Nerman’s theory of general branching processes counted by random characteristics. Another 
proof can be found in [5|, using only elementary probabilistic tools, relying on fluctuation analysis 
of the process. 

Our purpose in this work is to investigate the behaviour of the error in the aforementioned 
convergence. Many papers studied the second order behaviour of converging branching processes. 
Early works investigate the Galton-Watson case. In m and m, Heyde obtained rates of convergence 
and get central limit theorems in the case of supercritical Galton-Watson when the limit has finite 
variance. Later, in [T], Asmussen obtained the polynomial convergence rates in the general case. 
In our model, the particular case when the individuals never die (i.e. Py = 600 , implying that the 
population counting process is a Markovian Yule process) has already been studied. More precisely, 
Athreya showed in [3], for a Markovian branching process Z with appropriate conditions, and such 
that e~^^Zt converges to some random variable W a.s., that the error 

Zi - 

converges in distribution to some Gaussian random variable. 

In the case of general CMJ processes, there was no similar result although very recent work 
of Iksanov and Meiners m gives sufficient conditions for the error terms in the convergence of 
supercritical general branching processes to be o{t^) in a very general background (arbitrary birth 
point process). Although our model is more specific, we give slightly more precise results. Indeed, 
we give the exact rate of convergence, and characterized the limit. Moreover, we believe that 
our method could apply to other general branching processes counted by random characteristics, as 
soon as the birth point process is Poissonian. 

The first step of the method is to obtain informations on the moments of the error in the a.s. 
convergence of the process. Using the renewal structure of the tree and formulae on the expectation 
of a random integral, we are able to express the moments of the error in terms of the scale function of 
a Levy process. This process is known to be the contour process of the splitting tree as constructed 
in m- The asymptotic behaviours of the moments are then precisely studied thanks to the analysis 
of the ladder height process associated to a similar Levy process and to the Wiener-Hopf factoriza¬ 
tion. The second ingredient is a decomposition of the splitting tree into subtrees whose laws are 
characterized by the overshoots of the contour process over a hxed level. Finally, the error term can 
be decomposed as the sum of the error made in each subtrees. Our controls on the moments ensure 
that the error in each subtree decreases fast enough compared to the growth of the population (see 
Section 0] for details). 

The first section is devoted to the introduction of main tools used in this work. The first part 
recall the basic facts on splitting trees which are essentially borrowed from milSllTl El |5]. The 
second part recall some classical facts on renewal equations and the last part gives a useful Lemma 
on the expectation of a random integral. Section [3] is devoted to the statement of Theorem 13.21 which 
is a GLT for the population counting process W- Section [4] details the main lines of the method. 
Theorem 13.21 is finally proved in Section El 
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2 Splitting trees and preliminary results 


This section is devoted to the statement of results which are constantly used in the sequel. The first 
subsection presents the model and states results on splitting trees coming from dZl El El 1211 [5]. The 
second subsection recalls some well-known results on renewal equations. Finally, the last subsection 
is devoted to the statement and the proof of a lemma for the expectation of random integrals, which 
is constantly used in the sequel. 


2.1 Splitting trees 

In this paper, we study a model of population dynamics called a splitting tree. We consider a 
branching tree (see Figured]), where individuals live and reproduce independently from each other. 
Their lifetimes are i.i.d. following an arbitrary distribution Py. Given the lifetime of an individual, 
he gives birth to new individuals at Poisson times with positive constant rate b until his death 
independently from the other individuals. We also suppose that the population starts with a single 
individual called the root. 

The finite measure A := bFy is called the lifespan measure, and plays an important role in the 
study of the model. 



Figure 1: Graphical representation of a Splitting tree. The vertical axis represents the biological 
time and the horizontal axis has no biological meaning. The vertical lines represent the individuals, 
their lengths correspond to their lifetimes. The dashed lines denote the filiations between individuals. 

In [T7|, Lambert introduces a contour process Y, which codes for the splitting tree. Suppose we 
are given a tree T, seen as a subset of M x (Ua;>oN^) with some compatibility conditions (see [HD- 
On this object, Lambert constructs a Lebesgue measure A and a total order relation ^ which can 
be roughly summarized as follows: let x, y in T, the point of birth of the lineage of x during the 
lifetime of the root split the tree in two connected components, then y Y x y belong to the same 
component as x but is not an ancestor of x (see Figure Ej). 
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If we assume that A(T) is finite, then the application, 

(/?:¥—>• [0, A (T)], 

X A({y I y ^ xj), 

is a bijection. Moreover, in a graphical sens (see Figure [2|), (p(x) measures the length of the part of 
the tree which is above the lineage of x. The contour process is then defined, for all s, by. 



Figure 2: In gray, the set {y gT \ y < x} 

Y, := Hk (s)) , 

where 11® is the projection from M x (Ufc>oN^) to M. 

In a more graphical way, the contour process can be seen as the graph of an exploration process 
of the tree: it begins at the top of the root and decreases with slope —1 while running back along the 
life of the root until it meets a birth. The contour process then jumps at the top of the life interval 
of the child born at this time and continues its exploration as before. If the exploration process does 
not encounter a birth when exploring the life interval of an individual, it goes back to its parent and 
continues the exploration from the birth-date of the just left individual (see Figure [3|). It is then 
readily seen that the intersections of the contour process with the line of ordinate t are in one-to-one 
correspondence with the individuals in the tree alive at time t. 

In the case where A(T) is inhnite, one has to consider the truncations of the tree above fixed 
levels in order to have well-defined contours (see [17] for more details). In |17j . Lambert shows that 
the contour process s G of a splitting tree which has been pruned from every part above 

t (called truncated tree above t), has the law of a spectrally positive Levy process started at the 
lifespan V of the root, reflected below t and killed at 0, with Laplace exponent iIj given by 

V’(x) = X — f (l — A((ir), x G M+. (2.1) 

J (0,oo] 
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Figure 3: One-to-one correspondence between the tree and the graph of the contour represented by 
corresponding colours. 

In particular, the Laplace transform of Py can be expressed in terms of 'll:, 

[ e-^^Fvidv) = 1 + (2.2) 

Jr+ 0 

The largest root of V', denoted a, characterizes the way the population expend. In this paper, we 
only investigate the behavior of the population in the supercritical case, when a > 0. In particular, 
using the convexity of V’ (see [IS]), this is equivalent to ■0^(0-!-) < 0. Now, since 

ip'lx) = 1 — f xe~^'^ bFv{dv), Vx G M+, (2.3) 

J M-j_ 

one can see that the condition '0^(0-!-) < 0 is also equivalent to have 6E [F] > 1 which is a more 
usual supercritical condition. In the supercritical case, the population grows exponentially fast on 
the survival event with rate a. According to ()2.2p . one can also see that 

[ e-“"Pu(du) = 1 - X- (2-4) 

Jr+ 0 
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As said earlier, an important feature of the contour process is that the number of alive individuals 
at time t equals 

Card{yW = t I s G M+}. 

This set is the number of times the contour process hits t. This allows getting, thanks to the theory 
of Levy processes, the law of the unidimensional marginals of the process {Nt, t € K_|_). Indeed, let 
Tt (resp. To) be the hitting time of t (resp. of 0) by the contour process Now, for any positive 
integer k, the strong Markov property entails that 

P(Art = A; I AT* > 0) = EjPiAV = t | s > 0} = A: | r* < tq) } = P* t | s > 0} = - l) . 

Once again, the strong Markov property gives 

Pt = t I s > 0} = A: - l) = P* (rt < To)Pt = A | s > 0} = A; - 2 ) 

= Pt in < Pi = A I s > 0} = 0). 

Now, using fluctuation identities for spectrally positive Levy processes (see Theorem 8.1 in [IB] for 
the spectrally negative case), we have that 


Pi in < To) = 1 - 


1 

wwr 


where W is the scale function of the Levy process whose Laplace exponent is given by (12.11) . The 
function W is the unique increasing function whose Laplace transform is given by 

TcWit) = [ e-^^W{r)dr = ^, A > a, (2.5) 

i(0,oo) W[t) 

where a is the largest root of iJj. 

From the discussion above, we see that Nt is a geometric random variable conditionally on 
{Nt > 0}. More precisely, for a positive integer k, 




(2.6) 

In particular. 

E [Nt 1 iVi > 0] = 1F(A). 

(2.7) 

Moreover, it can 

be showed (see [H]), that 



ENt = W{t) - VF*Pv'(A), 

(2.8) 

and 

PW>o) = i-™*), 

(2.9) 

where 

W*Fv{t):=[ W{t-s)Fv{ds). 
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For the rest of this paper, unless otherwise stated, the notation refers to P (. | > 0) whereas 

Poo refers to the probability measure conditioned on the non-extinction event (which has positive 
probability in the supercritical case). 

Finally, we recall the asymptotic behaviour of the scale function W{t) which is widely used in 
the sequel. 

Lemma 2.1. Thm. 3.21]) There exist a positive constant 7 such that, 

e-‘^^'{a)W{t) -1 = 0 (e-^*) . 

From this Lemma and (12.9p . one can easily dednce that 

Q: 

P (NonEx) = lim F (Nt > 0) = —, (2-10) 

>00 0 

where NonEx refer to the non-extinction event. 

To end this section, let us recall the law of large number for Nt. 

Theorem 2.2. There exists a random variable £, such that 

g 

e-^^Nt ,,, , , a.s. and in . 

t=too ip'[a) 

Moreover, under Poo, £ is exponentially distributed with parameter one. 


2.2 A bit of renewal theory 

The purpose of this part is to recall some facts on renewal equations borrowed from [9]. Let h : M —)• M 
be a function bounded on finite intervals with support in R_|_ and F a probability measure on M+. 
The equation 

F{t)= [ F{t-s)Tids) + h{t), 

Jm.+ 

called a renewal equation, is known to admit a unique solution finite on bounded interval. 

Here, our interest is focused on the asymptotic behavior of F. We said that the function h is 
DRI (directly Riemann integrable) if for any <5 > 0, the quantities 

n 

s sup f(t) 

j _0 ts[<5*,<5(*-l-l)) 


and 


n 


S inf 

“ te[<Si,<5(i-l-l)) 


fit) 


converge as n goes to infinity respectively to some real numbers and 


and 


lim = lim r < 00 . 
< 5^0 ^ < 5^-0 •' 


In the sequel, we use the two following criteria for the DRI property: 
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Lemma 2.3. Let h a function as defined previously. If h satisfies one of the next two conditions, 
then h is DRI: 

1. h is non-negative decreasing and classically Riemann integrable on M+, 

2. h is cddldg and bounded by a DRI function. 

We can now state the next result, which is constantly used in the sequel. 

Theorem 2.4. Suppose that T is non-lattice, and h is DRI, then 

lim F{t) = 7 / h{s)ds, 


-1 


with 

7 := f / s Lids] 

\Jm+ 

if the above integral is finite, and zero otherwise. 

Remark 2.5. In particular, if we suppose that T is a measure with mass lower than 1, and that 
there exists a constant a > 0 such that 


/ 

Jr. 


e°T{dt) = 1, 


then, one can perform the change a measure 

f{dt) = e"*r(dt), 

in order to apply Theorem \2.4\ to a new renewal equation to obtain the asymptotic behavior of F. 


(See for details). This method is also used in the sequel. 

2.3 A lemma on the expectation of a random integral with respect to a Poisson 
random measure 

Lemma 2.6. Let f be a Poisson random measure on M+ with intensity 9\{da) where 9 is a positive 
real number and A the Lebesgue measure. Let also ^Xu\ u G be an i.i.d. sequence of non¬ 

negative cddldg random processes independent of f,. Let also Y be a random variable independent of 
f and from the family (Xu\ u G M+j . If f,u denotes f ([0, n]), then, for any t >0, 


i>l 


E [ X^^^hyyu f{du) = [ P(y > u)9EXudu, 
J[ 0 ,t] Jo 

where {Xu, u G M+) = u G . In addition, for any t < s, we have 


E 


'[o,t] 


^{dv) [ Xi^-hY>u f{du) 

Jl0,s] 


= / 9E[xl]E{Y>u) du 


[ [ 9^EXuEXuF {Y > u,Y > v) dudv. 
Jo Jo 


+ 


8 





Proof. Since the proof the two formulas lies on the same ideas, we only give the proof of the second 
equation. 

First of all, let / : —>• M+ be a positive measurable deterministic function. We recall that, for 

a Poisson random measure, the measures of two disjoint measurable sets are independent random 
variables. That is, for A, B in the Borel cr-field of M+, f{A n B^) is independent from f,{B), which 
leads to 

E [^{A)C{B)] = Ef{A)EC{B) + Var^(^ n B). 

Using the approximation of / by an increasing sequence of simple function, as in the construction 
of Lebesgue’s integral, it follows from the Fubini-Tonelli theorem and the monotone convergence 
theorem, that 

E f f{u,v)f,{du)f,{dv)= f 6f{u,u)du+ f f 6‘^f{u,v) dudv. 

i[o,t]x[o,s] Jo Jo Jo 

Since the desired relation only depends on the law of our random objects, we can assume without 

loss of generality that ^ is defined on a probability space (fl, T”, P) and the family (xi^\ s G E_|_) 

\ / 2>1 

is defined on an other probability space Then, using a slight abuse of notation, we define 

^ on n X Q by = fco, and similarly for the family X. 

Then, by Fubini-Tonneli Theorem, with the notation ([0, u]). 



But since the X^J are identically distributed and ^ is a simple measure (purely atomic with mass 
one for each atom) we deduce that, if u and v are two atoms of = Q if and only if u = u, 

which implies that 


f Xp\Q)Xp\Q)¥ (doj) 

Jn 


EXyEX^u, u ^ V, 
EXl, u = v, 


The result follows readily, and the case with the indicator function of Y is left to the reader. 


□ 


3 Statement of the theorem 


The a.s. convergence stated in Section 12.11 suggests to study the second order properties of this 
convergence to get central limit theorems. We recall that the Laplace distribution with zero mean 
and variance is the probability distribution whose characteristic function is given by 


A G M I-)- 


1 

1 + 
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It particular, it has a density given by 

™ 1 1^1 
X G M !->• —e . 

2 a 

We denote this law by £ (0 ,cj^). We also recall that, if G is a Gaussian random variable with 
zero mean and variance cr^ and E is an exponential random variable with parameter 1 independent 
of G, then y/SG is Laplace £ (O, cr^). 

Before stating the main result of the paper, let us recall the law of large number for W- 

Theorem 3.1. In the supercritical case, that is 6 E[I/] > 1, there exists a random variable £, such 

that ^ 

e~°‘^Nt , a.s. and in L^. 

t^oo ip'[a) 

In particular, under Poo; E is exponentially distributed with parameter one. 

In this work we prove the following theorem on the second order properties of the above conver¬ 
gence. 

Theorem 3.2. In the supercritical case, we have, under Poo, 

e-t* U'{a)Nt - e“‘£) ^ £ (O, 2 - V^'(a)) . 

The proof of this theorem is the subject of Section [6l Note that, according to (12.311 . we have 

2 — = 1 -|- f ve~°''“ bFvidv) > 0. 

J 


4 Strategy of proof 


Let {Gn)n>i be a sequence of geometric random variables with respective parameter and (W)j>i 
a family of i.i.d. random variables with zero mean independent of (Gn)„>^. It is easy to show 
that the characteristic function of 


Gn 


Zn ■— Wj, 

vn ^ 


2=1 


(4.1) 


is given by 

1 + + Onil) ’ 


(4.2) 


from which we deduce that converges in distribution to £(0,EW^). 

If we suppose that the population counting process N" is a Yule Markov process, it clearly follows 
from the branching property that, for s < t, 


Ns 




2 = 1 


(4.3) 
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where the family is an i.i.d. sequence of random variables distributed as Nts and inde¬ 
pendent of Ng. Moreover, since Ng is geometrically distributed with parameter taking the 

renormalized limit leads to, 

Ns 


lim =:£ = e”"" V Si, 

t —>-00 ^ 

i=l 

where Si,... ,Sn^ is an i.i.d. family of exponential random variables with parameter one, and inde¬ 
pendent of Ng. Hence, 


Ns 

Nt - e^^S = (^Ni_g - , 

i=l 

is a geometric sum of centered i.i.d. random variables. This remark and (|4.1I) suggest the desired 
CLT in the Yule case. 


Remark 4.1. Let N he a integer valued random variable. In the sequel we say that a random vector 
with random size form an i.i.d. family of random variables independent of N, if and only 

if 


where (Xi 


i>l 


{Xi,...,Xm) = [Xi,...,Xm) , 

is a sequence of i.i.d. random variables distributed as Xi independent of N. 


However, in the general case, we need to overcome some important difficulties. First of all, 
equation (14.31) is wrong in general. Nevertheless, a much weaker version of (14.31) can be obtained in 
the general case. To make this clear, if u < t are two positive real numbers, then the number of alive 
individuals at time t is the sum of the contributions of each subtrees T (Oj) induced by each alive 
individuals at time u (see Figure IH). Provided there are individuals alive at time u, we denote by 
the residual lifetimes (see Figured]) of the alive individuals at time u indexed using that 
the ith individual is the ith individual visited by the contour process. Hence, 


Nu 

Ni = YNum, ( 4 . 4 ) 

i=l 

where (A7_u denote the population counting processes of the subtrees T{Oi) induced by 

each individual. The notation refers to the fact that each subtree has the law of a standard split¬ 
ting tree with the only difference that the lifelength of the root is given by O*. More precisly, we 
define, for alH > 1 and o E M+, the population counting process of the splitting tree con¬ 

structed from the same random objects as the ith subtree of Figure dl where the life duration of 
the first individual is equal to o. Hence, from the independence properties between each individuals, 
(o), t>u,o> is a family of independent processes, independent of (Oi)i<j<jv„) 

[Nl_^{o),t > u) has the law of the population counting process of a splitting tree but where the 
lifespan of the ancestor is o. Note that the lifespans of the other individuals are still distributed as 
V. From the discussion above, it follows that the family of processes {Nl_^ (Oj), t > are 

dependent only through the residual lifetimes and the law of {Nt (Oi), t E M+) under 
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T(03) 


T(Oi) 


J 


T(04) 

O 4 


u 


Figure 4: Residual lifetimes with subtrees associated to living individuals at time u. 


is the law of standard population counting process of splitting tree where the lifespan of the root is 


distributed as Oi under P^. 

Unfortunately, the computation of (I4.2p does not apply to (14.41) . This issue is solved by the 
following lemma, whose proof is very similar to one of Proposition 5.5 of m- 

Lemma 4.2. Let u in M+, we denote by Oi for i an integer between 1 and Nu the residual lifetime 
of the ith individuals alive at time u. Then under P„, the family {Oi, i G [1, Al^]]) form a family 
of independent random variables, independent of Nu, and, expect Oi, having the same distribution, 
given by, for 2 < i < Nt, 



(4.5) 


Moreover, it follows that the family {Ns{Oi),s G is an independent family of process, i.i.d. 

for i >2, and independent of Nu- 


Proof. Let 


a family of independent Levy processes with Laplace exponent 


0<i<N, 


'tp{x) = X — f (1 ~ 6 ^{dr), X G M+, 

v/(0,oo] 

conditioned to hit {u, 00 ) before hitting 0, for i G {0,..., Nu — 1}, and conditioned to hit 0 first for 
i = Nu- We also assume that. 



and 


Y^^=u, iG IV4. 
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Figure 5: Reflected JCCP with overshoot over t. Independence is provided by the Markov property. 
Now, denote by r* the exit time of the iih. process out of (0, u) and 

n—1 

r„ = ^Ti, n G {0,..., + 1} . 

i=0 

Then, the process defined, for all s, by 

Nu 

Ys = Y,Ys-tM<s<T.^., 

i=0 

has the law of the contour process of a splitting tree cut above u. Moreover, the quantity Y^-^ — Y^-^- is 
the lifetime of the iih alive individual at time t. The family of residual lifetimes (Oi)i<j< 7 v„ then 
the same distribution as the sequence of the overshoots of the Y above u. Thus, the Markov property 
ensures us that (Oj, i G |2,Nu]]) is an i.i.d. family of random variables. The Markov property also 
ensures that Oi is independent of the other Oj’s. 

It remains to derive the law of Oj. Let T be a Levy process with Laplace exponent V’- We denote 
by T+ the time of first passage of —Y above u and Tq the time of first passage of —Y below 0. Then, 
for all i >2, 

Pn {Oi G dx) = Po {^-Y^- G dx I Tg" < . 

On the other hand. Theorem 8.7 of |16] gives for any measurable subsets A C [0,u], B C (0, —oo), 

P„ (-U- £ B. -U- _ £ a) = p_v (B - y) 
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The result follows easily from 


IP ('^0 < -^^) = 1 


1 


W{u)' 


□ 


Remark 4.3. It is important to note that the law of the residual lifetimes of the individuals considered 
above depends on the particular time u we choose to cut the tree. That is why, in the sequel, we may 
denote O) ' for Oi when we want to underline the dependence in time of the law of the residual 
lifetimes. 


In addition, as suggested by (j4.2p . we need to compute the expected quadratic error in the 
convergence of N^, 


E 





which implies to compute EA^^T. 

Although, this moment is easy to obtain in the Markovian case, the method does not extend 
easily to the general case. One idea is to characterize it as a solution of a renewal equation in the 
spirit of the theory of general CMJ processes. 

To make this, we use the renewal structure of a splitting tree: the splitting trees can be constructed 
(see m ) by grafting i.i.d. splitting tree on a branch (a tree with a single individual) of length Vg 
distributed as V. Therefore, there exists a family t G of i.i.d. population counting 

processes with the same law as {Nt, t € M+), and a Poisson random measure ^ on M+ with intensity 
b da such that 

Nt= N^tutv^yu C{du)+ tv^yt, a.s., (4.6) 

dm 


where ([0, u]). 

Another difficulty comes from the fact that unlike m, the quantities summed in (14.4p are 
time-dependent, which requires a careful analysis of the asymptotic behaviour of their moments. 

The calculus and the asymptotic analysis of these moments is made in Section 16.1.11 In Lemma 
16.11 we compute EiV^T, and then with Lemmas 16.21 and 16.41 we study the asymptotic behaviour of 
the error of order 2 and 3 respectively. Section 16.1.21 is devoted to the study of the same questions 
for the population counting processes of the subtrees described in Figured (when the lifetime of the 
root is not distributed as P). Finally, Section W?2\ is devoted to the proof of Theorem 13.21 

One of the difficulties in studying the behaviour of the moments is to get better estimates on the 
scale function W than those of Lemma EH This is the subject of the next section. 


5 Precise estimates on W using Levy processes 


Before stating and proving the result of this section, we need to recall some facts about Levy processes. 
We follow the presentation of m- First, we recall that the law of a spectrally positive Levy process 
{Yf, t G M+) is uniquely characterized by its Laplace exponent -0, 


V’y (A) = logE 



A G M. 
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which in our case take the form of (12.ip : 

V'y(A) = X- 


(l-e-™)6Py(dr), A G M+. 


(0,oo] 

In this section, we suppose that 10 = 0- For a such Levy process, 0 is irregular for (0, oo) and in this 
case the local time at the maximum [Lt, t G M) can be defined as 


nt 


Lt = '^e\ t 


G 




i=0 


where is a family of i.i.d. exponential random variables with parameter 1, and 

nt '■= Card{0 < s < t | Y* = supy^}, 

u<s 

is the number of times Y reaches its maximum up to time t. Finally, the ascending ladder process 
associated to Y is defined as 

Ht = sup Ys, t G M+, 


where G M+) is the right-inverse of L. It is known that LI is a subordinator whose values are 

the successive new maxima of Y. Conversely, in our case, the process (infs<tys, t G M_|_) is a local 
time at the minimum, denoted (^Lt, t G . The descending ladder process H is then defined from 
L as H was defined from L. 

We can now state, the celebrated Wiener-Hopf factorization which allows us to connect the 
characteristic exponent of Y with the characteristic exponents of the bivariate Levy processes 
{{Lt, Ht ), t G M+) and Ht^ , t G , respectively denoted by k and k. In our particular case, 

where Y is spectrally negative, we have 

f K{a, (3) = a, /? G M+, 

y k{a,f3) = (i)Y{oi) + 13, a,/3 gM+, 

where (py is the right-inverse of ipy. Taking a = 0 allows us to recover the Laplace exponent ipn of 
H from which we obtain the relation. 


V'y(A) = (A - (^y (0)) V’h(A). 

We have now all the notation to state and prove the main result of this section. 


(5.1) 


Proposition 5.1 (Behavior of W). There exists a positive non-increasing cadlag function F such 
that 


r^Oit 




t > 0, 


and 


lim e^^F{t) = 

>oo I g 


bvy-i < oo, 

otherwise. 
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Proof. Let Y'^ be a spectrally negative Levy process with Laplace exponent given by 

^«(A) = A - / (l- e-^^) e-^^b ¥v{dx). 

Jm.+ ^ ' 

It is known that has the law of the contour process of the supercritical splitting tree with lifespan 
measure Py conditioned to extinction (see [IZ!)- In this case the largest root of is zero, meaning 
that the process Y'^ does not go to infinity and that <^ytt(0) = 0. Elementary manipulations on 
Laplace transform show that the scale function W'^ of Y^ is related to W by 

W^t) = e-“*IL(t), t G M+. 

Let be the ascending ladder subordinator associated to the Levy process YK In the case 
where </>ytt(0) = 0, and in this case only, the scale function W'^ can be rewritten as (see m or use 
Laplace transform), 

W^{t) = [ P {hI < t) dx. (5.2) 

In other words, if we denote by U the potential measure of 

w^{t) = U[0,t]. 

Now, it is easily seen from (j5.ip that the Laplace exponent of takes the form, 

W = '*/’'(«) - [ (l - '^{dr), 

i[0,oo] ^ ^ 


where 


T(dr) = / e-‘^'’bFv{dv)dr = E [e-"^ly>^] bdr. 
J {r^oo) 


Moreover, 

T(M+) = 1-V^'(a), 

which mean that is a compound Poisson process with jump rate 1 — jump distribution 

J{dr) := dr, and killed at rate It is well known (or elementary by conditioning on 

the number of jumps at time x), that 

k>0 

Some calculations now lead to, 

U{dx) = ^T*^{dx). 

k>0 

From this point, since T is a sub-probability, U{x) := C/[0, x] satisfies the following defective renewal 
equation, 

U{x)= f U{x — u)T{du) + ls_^{x). 

J M-j_ 
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Finally, since 


e“^T(dx) = 1, 


and since, from Lemma 12.31 

t ^ U{t, oo), 

is clearly a directly Riemann integrable function as a positive decreasing integrable function. Hence, 
as suggested in Remark 12.51 

{U{R+) - U{x)) - 


1 


a/i 


with 


1 


^ = / re"'’T(dr) = - (6ER - 1), 

Jr+ “ 


if V is integrable. In the case where V is not integrable, the limit is 0. 

To end the proof, note using relation (j5.2p and the fact that H'^ is killed at rate tp^a) that, 

= 77 ^ -U{t,oo). 
ip'{a) 


□ 


6 Proof of Theorem 13.2 

We begin the proof of Theorem 13.21 by computing moments, and analysing their asymptotic be¬ 
haviours. A first part is devoted to the case of a splitting tree where the lifetime of the root is 
distributed as V whereas a second part study the case where the lifespan of the root is arbitrary (for 
instance, as the subtrees described by Figure 0]). 

6.1 Preliminary moments estimates 

This section is devoted to the calculus of the expectation of (W — We start with the simple 

case where the initial individual has life-length distributed as V. Secondly, we study the asymp¬ 
totic behavior of these moments. In Subsection 16.1.21 we prove similar result for arbitrary initial 
distributions. 

The expectations above are given with respect to P, however since W and T vanish on the 
extinction event, we can easily recover the results with respect to Pj by using (|2.10p and (|2.9I) (see 
Corollary 16.3p . 

6.1.1 Case 1/0 = 1/ 

We start with the computation of 
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Lemma 6.1 (Join moment of £ and Nt). The function t —>■ E [Nt£] is the unique solution bounded 
on finite intervals of the renewal equation, 


f{t) = f{t- M)6e-““P {V > u) du 


+ ahE[N.\i.(^j e-^^^{V > ■,¥ > v)dv^{t) 
+ a [ {V >t,V >v) dv, 


( 6 . 1 ) 


and its solution is given by 

(l + ^ - 6"“*) W{t) - (1 - e-“‘) Ws.Fv{t). 

Proof. As explained in Section HI 

At = / f,{du) + lv0>4, 

J[o,t] 

where f a Poisson point process with rate b on the real line, is a family of independent 

CMJ processes with the same law as N and V 0 is the lifespan of the root. Moreover, the three objects 
^ and Vn, are independent. 

It follows that, for s > t 


NtNs = [ Nt^N^iilv,>utv,>v f{du)adv) 

J[0,t]x[0,s] 

+ [ Tlt^u'^V0>u f.{du)lv^>s + I ^((iu)lyg>t + Iv0>tly0>s) 

J[0,t] J[0,d 

and, using Lemma [T6l 

EiVt = [ bE [Nt-uNs-u] P (P > u) du 
J[o,t] 

+ [ b'^E [Nt-u] E [Ns-v] E{V >u,V >v) du dv 

JlO,t]xlO,s] 

+ F{V > s) [ bE [Nt-u] du+ [ bE [A^,_J P (P > u, P > t) du + F{V > s). 

Then, thanks to the estimate IP(t) = O (e“*) (see Lemma 12.11 or IHTT]) and the convergence of 
IP(s)“^AtAs to Nt£ as s goes to infinity (since, by Theorem 12.21 converge in and using 
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Cauchy-Schwarz inequality), we can exchange limit and integrals to obtain, 


Ns 

lim EA^t-— 
s-J-oo w (s) 


= EN,£ = 




[ E [Nt-u£] e""“ P (1/ > u) bdu 

Jm 


+ 


X 


[0,t] X [0,oo) 


abE [Nt-u] {V > u,V > v) du dv 




+ [ ae-'^'’F{V > v,V >t) dv, 

J [0,oo] 


= -42{t) 


where we used that lim W{t) ^ENt = j. 


t—^OO 


Now, we need to solve the last equation to obtain the last part of the lemma. To do that, we 
compute the Laplace transform of each part of the equation. Note that, since W{t) = O (e"*), it is 
easy to see that the Laplace transform of each term of (16.ip is well-defined as soon as A > a (using 
Cauchy-Schwarz inequality for the first term). Now, using (12.21) . 


Ted 


'■G(A) = b [ {V >t)dt = b [ / Py (dv) dt 

j (t.oo) 


1 

A 


1 — e 


—Xv 


So, 

Then, 


TcG{X) = 1 - 


^ bPv (dv) = 1 - ^ 

ipiX + a) 


(t,oo) 

V'(A) 


( 6 . 2 ) 


A -)- Q; 


TcCiW = aTcEN,{X)Tc ( b / e’^^P iV > ■,¥ > v) dv ] {X) 

X 

V'(A) 


-l]Tc{a e-""P(l/>-,y>u)du (A). 


=.cc2(A) 


and, using (16.2p . we get 

TcC 2 {X) = a[ / e-^^P {V > t,V > v) dv dt = ]- ^ . 

Jr+ Jm.+ b \ X X j 

Finally, we obtain, 


TefiX) = TefiX) 1 - 


Hence, 


V>(A -F «) A / A 

A-I-a J VV’(A) 


- 1 


1 / 'ijj{X + a) _ V>(A) A 1 / il){X + a) _ V>(A) A 
b\ X X J b\ X X )' 


Tcf{X) = Y ^ 


1 


b y'tpW V'(A + a) 
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Finally, using (j2.5l) and 


bTc{Wi<Fv) (A) = 


(^(A)-6 + A) 

m ' 


allows to inverse the Laplace transform of / and get the result. □ 

Lemma [6TT] allows us to compute the expected quadratic error. 

Lemma 6.2 (Quadratic error in the convergence of Nt). Let £ the a.s. limit of ijj'{a)e~°‘^Nt. Then, 


Proof. Let 


lim e-“‘E (i>'(a)Nt - = ? (2 - . 

t^oo ^ ' b ^ 


[i := lim e“‘F(t), 

t—^OO 


where F is defined in Proposition 15.11 We have, using Proposition 15.11 and (12.41) . 


f e“‘ 

/ Wit-u)Fvidu) = -— 


1 - 


a 


[ e-^^Fv{du)+ [ (fi-e^^^-^^F{t-u))Fv{du) 

r (a) J{t,oc) J[o,t] ^ > 


6 ^^ / (y\ 

(i- )-p+„(i). 


e 


Hence, the expression of FNt£ given by Lemma 16.11 can be rewritten, thanks to Lemmas 15.11 as 

2ae"* a { 1 


FNt£ = 


+ h] + o(l), 


bfj'{a) b \'if'{a) 

Using ()2.6p and (|2.9p in conjunction with Proposition 15.11 we also have 

2a^ a 


e-^^FN'^ = 2, 


bip^aY b'il)'{a) b'ip'{a) 
Hence, it finally follows from (|6.3I) and ()6.4p that 


+ 0 ( 1 ). 


(6.3) 


(6.4) 
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at 


e-"*E {■4)'{a)Nt - e“*f) = f)'{af e-'^^FNf - 2i;'{a)FNt£ + 

= -2^V''(a) - _p 2^ (l + 'ip'{a)fi) + o(l) 

/y 


It is worth noting that, using (12.91) and the method above, we have the following result. 
Corollary 6.3. IFe have 

1 


□ 


’ {Nt >0) a 




a 


which leads to 


2e"* 1 

FtNt£ = . . - jrry - 3+ + o(l). 


(6.5) 

( 6 . 6 ) 
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Our last estimate is the boundedness of the third moments. 

Lemma 6.4 (Boundedness of the third moment). The third moment of the error is asymptotically 
bounded, that is 


E 




= 0 ( 1 ). 


Proof. We dehne for all t > 0, as the number of individuals alive at time t which have an inhnite 
descent. According to Proposition 6.1 of [5], N°° is a Yule process under P^o. 

We have 


E 


— i 
62 


< 


— f 
62 


3' 


N?° - 6"*T 

3" 


+ 8E 





6 2 



Now, we know according to the proof of Theorem 6.2 of [5] (and this is easy to prove using the 
decomposition of Figure 0]) that N°^ can be decomposed as 

Nt 


Nr = Y.Bf\ 


i=l 


where ( ) is a family of independent Bernoulli random variables, which is i.i.d. for i > 2, under 

V / i>l 
P^. Hence, 


E. 




62 ‘ 


< 6"i"*Ef 


Nt 


V i=l 




(t) 


Since, it is known from the proof of Theorem 6.2 of [5] that 

EBf = + O (e-“‘) , 

it is straightforward that 


E, 


fjfafNt - 


—f 
6 2 


is bounded. 

On the other hand, we know that a Yule process is a time-changed Poisson process (see for 
instance [2], Theorem III.11.2), that is, if Pt is a Poisson process independent of £ under Pqo, 


E 


^oo _ ^atg 

3' 

= Eoo 



3" 






6 2 




6 2 



P(NonEx). 


Now, using Hlder inequality, it remains to bound 

^ /V, . • . 

£(e“‘-l) 


En 


Pr(pat_T\ — 6“*T 


9Lf 

eP 


= 6"2“‘ / E„ 


(P3,(ect_i) - 


6 ^dx. 


Finally, for a Poissonian random variable X with parameter v, straightforward computations give 


that E 


{X — u) = -|- V, which allows us to end the proof. 


□ 
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6.1.2 Case with arbitrary initial distribution 

In order to study the behavior of the sub-splitting trees involved in the decomposition described 
in Figure 01 we investigate the behaviour of a splitting tree where the ancestor lifelength is not 
distributed as V, but follows an arbitrary distribution. Let S be a random variable in (0, oo], giving 
to the life-length of the ancestor and by N(E) the associated population counting process. 

Using the decomposition of N(E) over the lifespan of the ancestor, as described in Section 01 we 
have 

Nt(E) = / ^(du) + IsM, (6.7) 

J 

where is a family of i.i.d. CMJ processes with the same law as N independent of E and 

as described in section 01 Let, for all i>l,£i be 

£i := lim ip', a.s, (6.8) 

t—)-oo 

and, let £ (H) be the random variable dehned by 

£{E):= [ e(d^). (6.9) 

J [0,oo] 

Lemma 6.5 (First moment). The first moment is asymptotically hounded, that is 

W.{fi'{a)Nt{E) - e^^£{E)) = 0{l), 


uniformly with respect to the random variable H. 

Proof. Using Lemma ITUl (|6.7p and (|6.9I) with have 

E {fi'{a)Nt{E) - e“‘<S(“)) = J (fi'{a)ENt-u - ) e-““P {E > u) bdu, 


which leads using (|2.8p and (|2.10l) to 

E (fi'(a)Nt(E) - e^^£(E)) = [ (fi'(a)W(t - u) - fi'(a)W * Py(t - u) - e-““P (E > u) bdu. 

b 


— •If. — IJ 


( 6 . 10 ) 


We get using Proposition 15.11 and (|2.4p , 

/,=e“*-V^'(a)e“*F(s)-e“^(l-^) 

+ fi'{a) [ F{s — v)Fv{dv) + e' 

J [0,sl 


/ e-"^Py(du) - 

(s,oo) ^ 


=e 


[ e-^^Fv{dv)+o{l). 

J (s,oo) 


Hence, {Is)s>o is bounded. The result, now, follows from (I6.10p . 


□ 
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Lemma 6.6 (L^ convergence in the general case). ip'{a)e converge a.s. and in to £ (H), 

and 

lim e-“*E(V>'(a)A^i(H)-e“*£:(H))^ = ^ (2-V’'(a)) f e-^^¥ {E > s) bds, 

t^oo ' b ^ ^ 7r^ 

where the convergence is uniform with respect to H in (0, oo]. In the particular case when E follows 
the distribution of given by (j4.5p . we have, for 0 < /3 < ^, 

^lim = (2 - 

Proof. From (j6.7p and (j6.9p . we have 


(e-»V(a)JV,(E)-«(=))" = 


and, using Lemma [T6l 




£(u)) e f.idu) + e 

( 6 . 11 ) 


E{f;'{a)e-^^NtiE)-£{~)f 

+ e-2“*P (H > t) + 2e-“‘ElH>t [ U'ia)e-^d-u)^(du), 



— £{v)^ e "(“+'“)p (E > u,E > v) bdu dv 
e-““P (H > u, “ > t) bdu. 


Moreover, since, 

V^'(a)Ee-"*iVt -£ = 0 (e""*) , 

this leads, using Lemma [631 to 

lim e"*E (e-^^f;\a)Nt{E) - £{E)Y = %{2- [ e-"“P {E > u) bdu. 

t^oo ' b ^ 

Now, we have from (14.51) and Lemma l2.ll 

lim P„ (O 2 > s) = lim [ — ^P (V > s + y) bdy = [ e“"^P (V > s + y) bdy. 

u^co u^oo W{u)-1 

It follows then from Lebesgue theorem that, 

lim [ e~“^P/ 3 t {O 2 > s) bds = ^ ^ . □ 

t^oo _/r a 
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Lemma 6.7 (Boundedness in the general case.). The error of order 3 in asymptotically bounded, 
that is 

e-i"*E |V^'(a)A^i(H) - e^^£{E)f = O (1), 

uniformly w.r.t. H. 

Proof. Rewriting N{E) and £ (H) as in the proof of Lemma 16.61 we see that, 


e-i‘ E |V’'(a)iVt(E) - = e”!* E 

/ U'ia)Niii - ts>u adu) + ^P'ia)^E>t 

4 [04] ^ ^ 

< 8E 

[ y(a)y“j-e“ 

4[04] ^ 


+ 8V’'(a)e“^*P (H > t)^ 


We denote by I the first term of the r.h.s. of the last inequality, leading to 
3 

e~^''^lE>sif.{dsi)f,{ds2)f,{ds3) 


3 / tt \ 3 3 

g-i(i-Sj) (,/,// - N ArW,') 


< 8E / y 


< 


^ i=l 

24E J e-5(‘-“) (^^'(a)iVyj ^ e-^^ts>u^{du) e-^^^{du)^ 


< 24E 


'[04] 




with ^ 

li,{du) = ( f e“2®^(ds)| ^[du). 

\J[o,t] J 

Now, since /x is independent from the family and (^(i)), an easy adaptation of the proof of 

Lemma 1^61 leads to 

e“i‘ E |V’'(a)iVt(E) - <24E f E e-^^-u) (^^\a)Nt-u - 


404] 

+ 8'i/i'(Q;)e“2*P (E > t) 


e 2 “ 1 e:>« h{du) 


Using Lemma 16.41 to bound 


E 


(^Nt-u - e’^d-u)g'^ 


in the previous expression, finally leads to 


e-i* E \i;\a)Nt{E) - e^^£{E)f <C ^E , 


for some real positive constant C. 


□ 
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6.2 Proof of Theorem 13.21 


We fix a positive real number u. From this point, we recall the decomposition of the splitting tree 
as described in Section U] (see also Figure|l|). We also recall that, for all i in {1,..., Ny}, the process 
(IV] (O,), SGM+) is the population counting process of the (sub-)splitting tree T (Oi). 

As explained in Section 01 it follows from the construction of the splitting tree, that, for all i in 
{1,... ,Nu}, there exists an i.i.d. family of processes independent from Ny with the same 

law as {Nt, t € M_|_), and an i.i.d. family random measure independent from Ny and 

from the family with same law as such that 

Ni (O,) = / + to,>u Vt G M+, Vi G {1,... , Ny} . (6.12) 

J[o,t] 

As in (16.91) . we define, for all i in {1,..., Ny}, 

8 {Oi) := / f (6.13) 

where Sij := lim V’^(a)e““*iV]’A 

’ t^OO 

Hence, it follows from Lemma 16.61 that e~°‘^Nl (Oj) converges to E (Oj) in L^. 

Note also that, from Lemma 01^ the family (A"] (Oj), t G is i.i.d. and independent 

from Ny under P„, as well as the family [£ {Oi))2<i<Nu sense of Remark 14.11) . Note that 

the law under of the processes of the family (A] (Oj), t G is the law of standard 

population counting processes where the lifespan of the root is distributed as O 2 under P^ (except 
for the first one). 


Lemma 6.8 (Decomposition of £). We have the following decomposition of £, 

Nu 

f = e-““^fi(Oi), a-5- 

i=l 

Moreover, under Fy, the random variables {£i (Oi))j>i (defined by ()6.13p ) are independent, indepen¬ 
dent of Ny, and identically distributed for i >2. 


Proof. Step 1: Decomposition of £. 

For all t in M_|_, we denote by Nf^ the number of individuals alive at time t which have an infinite 
descent. For all i, we define, for all t > 0, N}f° (Oi) from T (Oj) as Nf° was defined from the whole 
tree. Now, it is easily seen that 

Ny 

Nr = ^Nffy{ 0 ,). 

i=l 

Hence, if (Oj) converges a.s. to £ (Oi), then 


lim e""*AA 

t—>-oo 


Ny 


Ny 


lim e-““ y (O*) = 6"““ V £ (Oi). 

t—^OO ^ ^ 


i=l 


2 = 1 
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So, it just remains to prove the a.s. convergence to get the desired result. 

Step 2: a.s. convergence of N°° (Oj) to £ {Oi). 

For this step, we fix i S {1,..., A^^}- 

In the same spirit as ()6.12p (see also Section^]), it follows from the construction of the splitting tree 
T (Oi), that there exists, an i.i.d. (and independent of Nu) sequence of processes s G M+ 

with the same law as (A^f°, t G M_|_) (under P), such that 

Nr m = [ io.>u + io,=oo, yt > o. 

Now, it follows from Theorem 6.2 of [5], that for all j. 


i>i 


lim e = £ij, a.s., 

t —^OO 

where £i j was defined in the beginning of this section. Let 


r := sup 


teK+ 


t ’ 


yj > 1 , 


and 


C := sup e-'^^Nr. 
tm+ 


Then, the family is i.i.d., since the processes are i.i.d, with the same law as C. 

Hence, 

/* r 

(6.14) 


J[0,t] J[0,t] 




It is easily seen that E [C] = P (NonEx) Eqo [C]. Now, since, from Proposition 6.1 of [5], Nr is a 
Yule process under P^o (and hence e~'^^Nr is a martingale), Doobs’s inequalities entails that the 
random variable C is integrable. Hence, the right hand side of the (|6.14ll is a.s. finite, and we can 
apply Lesbegue Theorem to get 


lim 


t—>oo 


m = [ 

dm 


£ Ai) e T{du) = £ (Oi ), a.s., 


where the right hand side of the last equality is just the definition of £ (Oi 


□ 


We have now all the tools needed to prove the central limit theorem for W- 


Proof of Theorem \3.^ Let u < t, two positive real numbers. From Lemma 16.81 and section [H we 
have 

Nu 

Nt = ^N^\m 


2=1 


and 


Nu 


^at£ = 


2=1 
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Then 


i^'{a)Nt - ^ gv (O,) - (0^) 


-t 
e 2^ 


i=l 


gf(t-«)gf« 


(6.15) 


Using Lemma IT^ we know that, under P„, (^Nl_^{Oi),t > are independent processes, i.i.d. 

for i > 2 and independent of N^- Let us denote by (p and (p the characteristic functions 


and 


(/j(A) := E 


^(A) := E 


exp iA 



i^'{a)Nl^ (02)-e“(*-“)T2(02)' 


Q ^{ t - u ) 


exp i\ 



(Oi)-e“(‘-“)Ti (Oi)' 


62 


{t-u) 


A G 


A G 


It follows from (j6.15p and Lemma 14.21 that. 


E„ 


exp iX— 


62 ' 










E„ 




e 2 


Nu 


Since Nu is geometric with parameter W(u) ^ under P, 

A 


E„ 


exp iX— 


—f 

62 






W{u)-^p 






Using Taylor formula for p, we obtain, 

E. 


, . U'(a)A^i-6“*T 
exp ^A— 


6 2 


= ^ 


X 


ef“y D{X,t,u) 


where, 

D{X,t,u) = W{u) 


- {W{u) - 1) ^l + iAE 


V^'(«)iV,T„(02)-e“(‘-“)T2(02) 




-^E 

2 


'i^'{a)NU (02)-e“(‘-“)T2(02)' 


e2A" 


-\- tj li) 




62 ' 


^P'ia)NU (O2) - e“(‘-“)T2 (O2) 


62 




2 e““ 


'^,^^^^{ 02)-£2 {O 2 ) 


ef 


- (lU(u) - l)i?(A,t,rt), 
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with, for all e > 0 and all A in (—e, e) 


\R{X, t, u)| < sup 

Ag(—e,e) 




'■/-'HiVt. (02)-e“l‘-“>f2(02)' 


^3 p 2 q 


< E I " ^ ^ J,. . -1 - - -< Ce^e-2“, 

(6.16) 

for some real positive constant C obtained using Lemma 16.71 

From this point, we set u = j5t with 0 < /3 < ^. It follows then from the Lemmas 16.61 and 14.21 
that 


lim E 

t—)-CxD 


•pt 


(O 2 ) - e“(*-/5b£:2 (O 2 )' 


efh-/3d 


='tjj'(a) (2 -'tjj'(a)) . (6.17) 


Moreover, we have from Lemma 16.51 and since /3 < ^, 

lim VF(/3f)e-t*E \Tp'{a)Ni (O 2 ) - e“*f 2 ( 02)1 = 0. 

t —>-00 

Finally, the relations (|6.16p . (16.171) and (|6.18l) lead to 

1 


(6.18) 


lim Ept 

t —>00 


exp iX 


Nt - 

27 

62 


l + ¥(2-V’'(a)) 


To conclude, note that, 


E 


pt 


exp iX 




62' 


-E,, 


exp ( iX 


Nt - 6“*T 


62 ' 


E 


i\ 




^Npt>0 


LnohEx \ 


P {Njst >0) P (NonEx) / 


< E 


^Npt>0 

1 NonEx 

_ P {X^pt > 0) 

P (NonEx) 


goes to 0 as t goes to infinity. This ends the proof of Theorem 13.21 


□ 
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